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Abstract 

In this paper, we study the perfect and the insulated conductivity problems with multiple 
inclusions imbedded in a bounded domain in R n ,n > 2. For these two extreme cases of the 
conductivity problems, the gradients of their solutions may blow up as two inclusions approach 
each other. We establish the gradient estimates for the perfect conductivity problems and an upper 
bound of the gradients for the insulated conductivity problems in terms of the distances between 
any two closely spaced inclusions. 

Introduction 

In this paper, a continuation of [5], we establish gradient estimates for the perfect conductivity problems 
in the presence of multiple closely spaced inclusions in a bounded domain in M" (n > 2). We also 
establish an upper bound of the gradients for the insulated conductivity problems. For these two 
extreme cases of the conductivity problems, the electric field, which is represented by the gradient of 
the solutions, may blow up as the inclusions approach to each other, the blow-up rates of the electric field 
have been studied in [TJ [3J [51 1191 120] . In particular, when there are only two strictly convex inclusions, 
and let e be the distance between the two inclusions, then for the perfect conductivity problem, the 
optimal blow-up rates for the gradients, as e approaches to zero, were established to be e -1 / 2 , (e\ lne|) _1 
and e _1 for n = 2, 3 and n > 4 respectively. A criteria, in terms of a functional of boundary data, for 
the situation where blow-up rate is realized was also given. See e.g. the introductions of [5] and for 
a more detailed description of these results. More recently, Lim and Yun in [T5] have obtained further 
estimates with explicit dependence of the blow-up rates on the size of the inclusions for the perfect 
conductivity problem (see also [1] for results of this type), and H. Ammari, H. Kang, H. Lee, M. Lim 
and H. Zribi in [2] have given more refined estimates of the gradient of solutions. 

The partial differential equations for the conductivity problems arise also in the study of composite 
materials. In R 2 , as explained in [14] , if we use the bounded domain to represent the cross-section of a 
fiber-reinforced composite and use the inclusions to represent the cross-sections of the embedded fibers, 
then by a standard anti-plane shear model, the conductivity equations can be derived, in which the 
electric potential corresponds to the out-of-plane elastic displacement and the electric field corresponds 
to the stress tensor. Therefore, the gradient estimates for the conductivity problems provide valuable 
information about the stress intensity inside the composite materials. 

When conductivities of the inclusions are away from zero and infinity, the boundedness of the gra- 
dients were observed numerically by Babuska, Anderson, Smith and Levin [4]. Bonnetier and Vogelius 
[6] proved it when the inclusions are two touching balls in M 2 . General results were established by Li 
and Vogelius |14j for second order divergence form elliptic equations with piecewise smooth coefficients, 
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and then by Li and Nirenberg [13] for second order divergence form elliptic systems, including linear 
system of elasticity, with piecewise smooth coefficients. See also [H) and [H] for related studies. 
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1 Mathematical set-up and the main results 

Let SI be a domain in M. n with C 2 ' a boundary n > 2, < a < 1. Let {Di} (1 < i < m) be m strictly 
convex open subsets in SI with C 2:Q boundaries, m > 2, satisfying 



the principal curvature of dDi > Kq, 
Sij :=dist(A,-Di) >Q, (i^j) 

dist(Di,50) > r , diamfQ) < — , 

TO 



(1.1) 



where «o ; ro > are universal constants independent of We also assume that the C 2a norms of 

dDi are bounded by some universal constant independent of {£ij}. This implies that each Di contains 
a ball of radius Tq for some universal constant r$ > independent of {sij}. 

We state more precisely what it means by saying that the boundary of a domain, say SI, is C 2a for 
< a < 1: In a neighborhood of every point of dSl, dSl is the graph of some C 2 ' a function of n — 1 
variables. We define the C 2 a norm of dfl, denoted by ||9f2||c2,c<, as the smallest positive number ^ 
such that in the 2a— neighborhood of every point of 9f2, identified as after a possible translation and 
rotation of the coordinates so that x n = is the tangent to dft at 0, dil is given by the graph of a C 2,a 
function, denoted as /, which is defined as \x'\ < a, the a— neighborhood of in the tangent plane. 
Moreover, ||/||c 2 -=(|x'|<a) < 

Denote 

Q := ft\U™iA- 

Given ip G C 1,Q (9J7), the conductivity problem can be modeled by the following equation: 



div(ak{x)Wuk) = in O, 
Uk = <p on dfl, 



(1.2) 



where k = (k\, . . . , k m ) and 



a k {x) 



ki e (0, oo) in Di, 
1 in f2. 



(1.3) 



The existence and uniqueness of solutions to the above equation is well known. Moreover, we have 
||wfe||/fi(o) < C\\ip\\ci-"(dfi) for some constant C independent of k. Therefore, by passing to a subse- 
quence, we have Uk — v "oo in H (SI) as fc, — > oo for all 1 < i < m, where itoo € H 1 ^) is the solution 
to the following perfect conductivity problem, 



( Au = 

u|+ = u| 
Vu = 

dDi dv 
t U = if 



where 



du 

dv 



= 



:= lim 

+ t->o+ 



in 17, 

on dDi, (i = 1, 2, . . . , m), 
in Dj (i = 1, 2, . . . , m), 

(i = 1,2,.. . ,m), 

on 9f2, 

tt(a; + fi/) — u(x) 



(1.4) 
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Here and throughout this paper v is the outward unit normal to the domain and the subscript ± 
indicates the limit from outside and inside the domain, respectively. For the derivation of the above 
equation, readers can refer to the Appendix of [5], Note that the proof there is for ki = = • • • = k m , 
but it works also for the general case with modification. 

Since the high stress concentration only occurs in the narrow regions between the fibers, we only need 
to focus on those narrow regions. 
For i ^ j, denote 



and 



dist(x y ,Xy) = dist(A,£>j) = £%j > 0, x y G 8D i: Xy € dD h 



c ij -~ 2^ X ^ Xi i^' 



It is easy to see that there exists some positive constant 5 < 4 which depends only on kq, tq and 
{||9£)i||c!3.o>}) but is independent of {ey} such that 



if Eij < 26, B(x®j,25) only intersects with Di and Dj. 



(1.5) 



Denote 



1 

1 



e| lne| 
1 

£ 

Then we have the following gradient estimates for the perfect conductivity problem 



for n = 2, 
/or n = 3, 
/or ?i > 4. 



(1.6) 



Theorem 1.1 Let ft, {Di} C K™, n > 2, {ey} &e defined as in ip G L°°(dft), 6 be the universal 

constant satisfying (| 1 . 5\ . Suppose u M € H 1 (r2) ?s i/ie solution to equation then for any £y < <5. 

we have 

l|Vuoo|| L! »(Q nBCx o. i5)) < Cpn(£»j)IMU°°(8n) 

where C is a constant depending only on n, Kq, Tq, {\\dDi\\c^,^} , but independent of Eij. 

Note that if £y > 5, by the maximum principle and the boundary estimates of harmonic functions, 
we immediately get ||VKc»||£oo(£2 nS (, £ .o s \\ < CIMI^^-, < C||y||z,°°(ar2) • Here we have used the fact 
that is constant on each dDi. Then by Theorem 11.11 and standard boundary Schauder estimates, 
see e.g. Theorem 8.33 in [S], we have the global gradient estimates of in Ct. 

Corollary 1.1 Let O, {Di} C M. n , n > 2, {s^} be defined as in e := min£y > 0, and tp G 

C 1:Q (9il), < a < 1, and let G H (O) be the solution to equation jl.4\ )- Then 

\\Vuoo\\ Loa( n) ^ c 7>n(e)IMIc?i.<»(an)- 

where C is a constant depending only on n, m, kq, ro, \\dCl\\c^,a, {\\dDi\\c^-"}, but independent of e. 

Remark 1.1 The proof of Theorem \1.1\ does not need Di and Dj to be strictly convex, the strict 
convexity is only used in a fixed neighborhood o/xy (The size of the neighborhood is independent of 
{eij}). In fact, our proofs of Theorem \T7l\ also apply, with minor modification, to more general situations 
where two closely spaced inclusions, Di and Dj , are not necessarily convex near points on the boundaries 
where minimal distance s is realized; see discussions after the proof of Theorem \l.l\ in Section 2. 

Next, we study the insulated conductivity problem. Similar to the perfect conductivity problem, 
the solution to the insulated conductivity problem is also the weak limit of Uk in H (O) as k approaches 
to 0. Here we consider the insulated conductivity problem with anisotropic conductivity. 
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Let Fl,Di C R", £ij be defined as in ip G C 1,a (dft), suppose A(x) := (a z:, (x)) is a symmetric 

matrix function in f2, where a u (x) G C Q (il) and for some constants A > A > 0, 

||a y |l CQ(I) <A, a y '(x)^>A|£| 2 , V£eR",xea 

Then the anisotropic insulated conductivity problem can be described by the following equation, 

' di(a ij dju) = in n, 

< tf'djUVi = on 9A(i = 1, 2, . . . ,m), (1.7) 

t u = ip on dtt. 

The existence and uniqueness of solutions to equation (|1.7J) are elementary, see the Appendix. 

As mentioned before, the blow-up can only occur in the narrow regions between two closely spaced 
inclusions. Therefore, we only derive gradient estimates for the solution to (|1.7p in those regions. 
Without loss of generality, we consider the insulated conductivity problem in the narrow region between 
D\ and D<i. Assume 

£ = dist(L>i,D 2 ) 
After a possible translation and rotation, we may assume 

(e/2,0') e 8D U (-e/2,0') e dD 2 . 

Here and throughout this paper by writing x = (xi, x'), we mean x' is the last n — 1 coordinates of x. 
We denote the narrow region between D\ and £>2 and its boundary on dD\ and dD 2 as follows 

0(r) :=f2n{xGR"||x'| < r} 

r+ :=&Din{xeR"||x'| <r} (1.8) 
r_ :=<9L> 2 n{xGR"||x'| < r} 

where r is some universal constant depending only on {||i9Z3i|| C 2,c t }. 
With the above notations, we consider the following problem, 



(1.9) 



J di{a ij dju) = in 0(r), 

\a ij djUVi = Q onr+ur_. 

Then we have: 

Theorem 1.2 If uq G i/ 1 (C(r)) is a weak solution of U.9]) . then 

|V.o(x)| < C ^ h ~™\ for all x G 0(1). (1.10) 

where C is a constant depending only on n, k q , r Q , A, A ; r and \\dDi\\ca,*(i = 1,2), but independent of 
e. 



Remark 1.2 Theorem also remains true for general second order elliptic systems, its proof is 
essentially the same as for the equations. 

A consequence of Theorem 11.21 is the following global gradient estimates for the insulated conductivity 
problem. 

Corollary 1.2 Let Q, {D t } C R n , {e. y } be defined as in / TO)) , e := miney > 0, and ip G C^ a (dO,), let 
uq G i/ 1 (r2) be the weak solution to equation \1. 7| ), then 

C 

||Vwo|| L =o ( o) < -j=\\ip\\ C i,<*(m)- (1-H) 
where C is a constant depending only on n, kq, tq, ||9Q||c2,a, {HS-DiHc 2 ."}; but independent of e. 
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Note that throughout this paper we often use C to denote different constants, but all these constants 
are independent of e. 

The paper is organized as follows. In Section 2 we consider the perfect conductivity problem and 
prove Theorem ll.il In Section 3 we show Theorem II .21 for the insulated case. Finally in the Appendix 
we present some elementary results for the insulated conductivity problem. 



2 The perfect conductivity problem with multiple inclusions 



In this section, we consider the perfect conductivity problem (jl.4[> . Note that from equation (|1.4p . we 
know that u = d on Di 1 1 < i < m, where {Cj} are some unknown constants. In order to prove 
Theorem 11.11 we first estimate \Ci — Cj\ for 1 < i ^ j < m, which later will allow us to control the 
gradient of u in the narrow region between Di and Dj . 

2.1 A Matrix Result 

To estimate \Ci — Cj\, the following proposition plays a crucial role. 

Let m be a positive integer, P = (pij) an m x m real symmetric matrix satisfying, 

(Al) Pij =pji<0 (i^j); 

m 

(A2) 0<n <p, :=^2pij <r 2 , 
i=i 

where r\ and r 2 are some positive constants. 

Remark 2.1 An m x m matrix P satisfying \pa\ > ^2 \Pij\ * s called a diagonally dominant matrix. 
Such a matrix is nonsingular, see (^41) and (A2) imply that the matrix P is diagonally dominant. 

Proposition 2.1 Let P = (pij) be an m x m real symmetric matrix satisfying (Al) and (A2), m > 1. 
For [3 € R m , let a be the solution of 

Pa = /3, (2.1) 



then 



aj \<m(m-l)^ f\ , (2.2) 

r l \Pij I + r l 



where \[3\ = max|/3i|. 

i 

Before proving the proposition, we introduce the following lemmas. 
Denote 

1(1) = {all / x I diagonal matrices whose diagonal entries are 1 or —1}, 
T e (l) = {I G I{1) \l has even numbers of —1 in its diagonal}, 
T {1) = {I S 1(1) \l has odd numbers of —1 in its diagonal}. 

Lemma 2.1 For any i£l and any I x I matrix A, I > 1, 

det(xI+lA) = 2 l - x (x l +detA); 

det(xI + lA) = 2 l -\x l -detA). 
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Proof: We prove it by induction. The above identities can be easily checked for 1 = 1. Suppose that 
the above identities stand for I = k — 1 > 1, we will prove them for I = k. Observe that the above 
identities hold when x = 0. To prove them for all x, it suffices to show that the derivatives with respect 
to x in both sides of the identities coincide. Since for any / S T(k), 

k 

(det (xl + 1 A))' = det ( xI + Ii A i) 

i=i 

where and Jj are the submatrices obtained by eliminating the ith row and the ith column of A and 
I respectively. 

Notice that if I runs through all the elements of l e (k), It will run through all the elements of I(k — 1) 
for every fixed i € {1, 2, . . . , k}, so we have 

(det(xI + IA))' 

7ez„ (fc) 
k 

= J2( det (xl + 7Ai) + Y det (xI + IAi)) 

i=1 7ex e (fc-i) 7ex (fc-i) 
fc 

= ^ (2 fc - 2 (a; fc - 1 + det A) + 2 k - 2 (x k - 1 ~ det A,)) (By induction) 

i=l 

= k2 k ~ 1 x k - 1 = 2 k - 1 (x k +detA)'. 

Therefore, we have proved the first identity. The second one follows from the first one by changing the 

sign of one row of A. 

As a consequence of Lemma 12. 1[ we have 

Corollary 2.1 Let A be an I x I matrix, if det (/ + I A) > for any I G 1(1), then \ det A\ < 1. 

Lemma 2.2 Given integers m > I > 1, let Q — (qij) be an m x I real matrix which satisfies, for 
J = 1,2,..., I, 

Let A be the set of all I x I submatrices of the above matrix Q and Si € A the matrix obtained from the 
first I rows of Q, then we have 

det Si = max I det SI. 

seA 1 1 

Proof: For any S £ A, by rearranging the order of its rows we do not change | detS"!. Thus we can 
treat S as a matrix obtained by replacing some rows of Si by some other rows of Q. Note that S and 
Si could have no rows in common, which means S is obtained by replacing all the rows of Si by some 
other rows of Q. 
Given any / € 2^(0 > we claim: 

det (Si +1S)>0 
Proof of the claim: There are two cases between Si and S : 

Case 1. Si and S have no rows in common. Then by (|2.3j) . we know that Si+LS is diagonally dominant, 
therefore det (Si + IS) > 0. 

Case 2. Si and S have some common rows, denote the order of these rows byl < ii < ••• < i s <l,l< 
s < I. If row i S[) of LS is opposite to row i SQ of S for some 1 < sq < s, then row i So of Si + LS is 0, 
therefore det (Si + IS) = 0. Otherwise row it of IS is the same as that of S and Si for any 1 < t < s, 
then we take out the common factors 2 in these rows when we compute det (Si + IS), thus we have 

det (Si + IS) = 2 s det (Si + IS), 



G 



where S is the matrix obtained by replacing row it of S by for any 1 < t < s. We know that Si + IS 
is diagonally dominant according to (|2.3p , then det (Si + IS) > 0, it yields that det (Si +IS) > 0. 
Therefore, the claim is proved. 
Since det Si > and 

det (Si + IS) = det (I + ISSf 1 ) det Si 
we have, by the claim, that for any I £ 1(1), 

det (7 + TSSf 1 ) > 

By Corollary |2.1i we have 

| det (SSf/ 1 )! < 1 

therefore 

det Si > | det SI. 
Now we are ready to prove Proposition 1 2. II 



Proof of Provosition \2. li For m — 1 the inequality is automatically true. For m — 2, we have, by 
Cramer's rule, 



Pi 


Pl2 




Pll 


Pi 




01 


Pi 




P22 




P21 


02 




02 


P2 


Pll 


Pl2 




Pll 


Pl2 




Pll 


Pl2 


P21 


P22 




P21 


P22 




P21 


P22 



Since ri < pi < r 2 by Condition (A2), 



01 Pi 

02 P2 



0iP2 - 02Pi < 2r 2 



On the other hand, by Condition (Al) and (A2) 

P1P22 -P2P12 >PiP22 > ri(n + \piz\). 



Pll 


P12 




Pi 


Pl2 


P21 


P22 




P2 


P22 



Therefore, Proposition 12. II for m = 2 follows from the above. 

For m > 3, we only estimate \ati — a 2 \ since the other estimates can be obtained by switching 
columns of P. 

Since a satisfies (|2.ip . by Cramer's rule, we have: 



ai — a 2 



01 


P12 




Pirn 






Pll 


01 ■■ 


Plm 


02 


P22 




P2m 






P21 


02 ■■ 


P2m 




Pm2 




Pmm 






Pml 




Pmm 


pn 


Pl2 




Plm 




Pll 


P12 ■ 


■ Plm 


P21 


P22 




P2m 




P21 


P22 ■ 


■ P2m 


Pml 


Pm2 




Pmm 




Pml Pm2 


Pmm 


01 


Pll ^ 


-Pl2 


Pl3 






Plm 






02 


P21 ^ 


'P22 


P23 






P2m 






03 


P31 ^ 


~P32 


P33 






P3m 






0m 


Pml ^ 


-Pm2 


Pm3 




Pmm 








Pll 


Pl2 




Plm 












P21 


P22 




P2m 












Pml 


Pm2 




Pmm 
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By adding the last (m — 2) columns of the matrix in the numerator to its second column, we have 



Pi Pi 


Pl3 


Pis 




Pi P2 


P23 


P2s 




P3 P3 


P33 


P3s 




Pm Pm 


Pm3 


Pm.rn 


detP 




Pll 


Pl2 ■ ■ 


Plm 


detP 




P21 


P22 ■ ■ 


P2m 






Pml 


Pm2 ■ ■ 


Pram 





Next we estimate the determinants of the above two matrices separately. 
Expanding det P with respect to the first column, we have 

m 

detP = ^ Pjl P jl 

3 = 1 

where Pji is the cofactor of pj\. 

Applying Lemma l2.2l to the m x (m — 1) matrix obtained by eliminating the first column of P, we know 
that, among the cofactors Pji, Pn > has the largest absolute value. Since pji — pij < (j ^ 1) and 
Pn > by condition (Al) and (^42), we have 



detP>J2Pji p n = Pi Pu- 

3 = 1 



For the same reason, we have 



Pn = 



P22 



P2r, 



Pm2 ■ ■■ P,i 



i=2 



P33 



Pm3 



P3r. 



Combining the above two inequalities and using condition {A\) and (A2), we have 



detP > pi y^P2j 

3=2 



P33 



P3r, 



Pl(P2 ~P2l) 



Pm3 ' ' ' Pram 
P33 ■■■ P3 



P33 



P3r, 



Pm3 ■ - Pr, 



> riQpul + ri) 



Pm3 



(2.4) 



By Laplace expansion, see e.g. page 130 of [T7], we can expand detP with respect to the first two 
columns of P, namely, 

Pii Pn 
Pi 2 P12 



det P = 



P 



ilia 12: 



(2.5) 



where 1 < ii < ii < m and P; 1 i 2 i2 is the cofactor of the 2nd-order minor in row and column 1, 2 
of P. 

Applying Lemma 12.21 to the m x [m — 2) matrix obtained by eliminating the first 2 columns of P, we 
know that, among all those cofactors, 



P33 



P3r, 



Pm3 - ■ Pr, 
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has the largest absolute value. Since < pi < r 2 by condition {A2), 

f f <2r 2 |/3|, 

then by (|2.5p . we have 

I det P| < m(m — l)r 2 
By dHJ) and dM]), we have 



P33 



Prn3 



P3r. 



Pr, 



\0i\ - a 2 \ 



IdetPl r 2 
< to(to — 1) 



I det PI 



2.2 Proof of Theorem CCD 

As in [5], we decompose Uoc into to + 1 parts: 



Woo = «0 + y^^CjVj, 



(2.6) 



(2.7) 



where Uj € P x (f2) (i = 0, 1, 2, . . . , m) are determined by the following equations: 
for i = 0, 

' Au = 
«o = 
u = 



in fi, 

on dDi , dD 2 , ■ ■ ■ dD Tl 
on dfl. 



for i — 1,2,..., to, 



f Avi = 

t'i = 1 
Vi = 
v «i = 



on dDi, 

on dDj, for j i, 
on dfl. 



(2.8) 



(2.9) 



Since Uoc satisfies the integral conditions in equation (|1.4|) . using the decomposition formula ()2.7j) . we 
know that the vector (C±, C 2 , . . . , C m ) satisfies the following system of linear equations 



/ an 012 



«2r 



where 



\ 


c 2 




b 2 


/ 








/ = 


1,2,. 


. . , nr 


), 



(2.10) 



(i = l,2,...,m). 



(2.11) 
(2.12) 



Similar to the two inclusions case in [3], we first investigate the properties of Vi (i = 0, 1, • • • , to), 
the matrix A = (fly) and the vector 6 defined by (|2. 1 1|) and (|2.12|) . Here we state the following lemma, 
for its proof, readers may refer to Lemma 2.4 in [5]. 

Lemma 2.3 For 1 < i, j < to, let a,ij and 6, be defined by i2.ll}) and h2.12\) , then they satisfy the 
following: 
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(1) an < 0, dij = a,ji > (i^j), 

(2) -C< o*i<-^> 

l<j<m 

(3) \bi\<C\\ V \\ L oo {dn)> 

where C > is a universal constant depending only on n, Kq, Vq, ||e?fi||cra,a, but independent of Eij. 

Remark 2.2 From property (1) and (2) in Lemma \2.3[ we know that A is diagonally dominant, there- 
fore it is nonsingular. 

Lemma 2.4 Let Vo,Vi(i — l,...,m) be the solutions of equations \2.8\) and V2.9f) respectively, 8 is 
the constant satisfying il. 5\) , then there exists a universal constant C depending only on n, m, Tq, Kq> 
\\dDi\\c2,a and \\dCl\\c*><* i but independent of {Eij} such that, 

(1) ||V«o|| ioo(fi) < C; 

(2) llV^H^^^.^na) < % < s i 

(3) |Vui| < C on n \ ( U jyi>ey <5 B{x% , 6)) . 

Proof: The proof of (1) is the same as the proof of Lemma 2.3 in [5]. Since HvjII^oo/q) = 1, 8 is the 
constant satisfying (|1.5p , then if < 8, then by (|1.5[) . we know that B{x^-, 8) only intersects with £), 
and Dj, and B(x^,S) is at least 6 away from other inclusions. Then (2) just follows from the maximum 
principle and standard boundary estimates for harmonic functions. For the same reason, to prove (3), 
we only need to prove ||Vui||jr )00 / B / a .o s)nh) < C li k,l ^ i and eui < 8. Without loss of generality, we 
assume k = 1, 1 = 2, i — 3. Let V3 be the solution of the following equation, 



'Av 3 = in fl\DiUD 3 , 

U3 = on dDi, 

«3 = 1 on dD 3 , 

W3 = on <9f2. 

Then we have V3 > v 3 on <9f2, by the maximum principle, U3 > v 3 in SI. Since V3 — v 3 = on dD\, we 
have _ 

dv dv 

But |V«3| < C on dDi n B(x\ 2 i $) by the boundary estimates of harmonic functions, then we have 

l|V«3||i~(a£) 1 nB(x? 3 ,a)) = ll-g^-IU^caDinsCx^^)) < ( 2 -13) 

Similarly, we have 

l|V«3|U«>(a£» 2 nB(a!; 3 ,i)) = H"g^"lli< : »(aD2nB(x; 2 ,5)) < ( 2 -i4) 

Furthermore, by gradient estimates and boundary estimates of harmonic functions, we have 

Since VU3 is still harmonic function on B(xi 2 ,8) fl O, by (|2.13[) . (|2.14[) and (|2.15p and the maximum 
principle, we have 

ll V ^lli-(nnB( a ? 9 ,«)) < C. 

Next, we derive some further estimates of A = ((%■). 
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Lemma 2.5 Let be defined as in V2.ll)) , then there exists a universal constant C > 0, depending 
only on n, Tq, Ko, \\dDi\\ci,a and ||90||c2,c, but independent of {£ij}, such that for 1 < i ^ j < m, 

C 1 1 C 

< an < -——==^=, < dij < ——, for n = 2, 



min£j fc C /nime ik ' CJe^ 



-C\ ln(mine ifc )| < a it < — — | ln(minej fc )|, — | Ihe^I < a i3 < C\ lney|, for n = 3, 
-C < an < — — , — < ay < C, for n > 4. 



l :i 

Proof: Without loss of generality, we assume i = 1, j = 2. The proof of the estimates for an is the 
same as that in Lemma 2.5, Lemma 2.6, and Lemma 2.7 in [5]. Here we prove the estimate for a\2- In 
the following, we use C to denote some universal constant depending only on n, ro, koj ||<9-Di|lc 2 ' Q an< i 
||<9£1|| C 2, Q , but independent of {e^ }. 

Notice that if £12 is larger than some universal constant, then the proof is trivial. Therefore, we 
can assume £12 < 8, where S < 1/4 is the universal constant satisfying (|1.5[) . By (|1.5p , we know that 
B(xi2,S) only intersects with D\ and D 2 . 

Denote 

:= dDi n B(x° 12 , S) (i = 1, 2), T 3 := dB(x° 12 , 6) \ (D 1 U D 2 ) 
Since B{x\ 2: 2<5) does not intersect with Di(i > 3) or dfl by (|1.5p . then 

dist(T 3 ,UZ 3 dDi) > 5, dist(T 3 ,dQ) > 5, 
by standard gradient estimates and boundary estimates for harmonic functions, we have 

||Vwi|U-(r,)<C (2.16) 

By Lemma l2~4l we have ||Vui||i,ao(a.D 3 \r 3 ) < C. 
Therefore, we have 

f dvi f dvi f dvi f dvi 
oia=/ -5-=/ -H-+/ -5-=/ «^ + 01 - 2.17 

By the harmonicity of v\ on B(xi 2 , S) D and (|2.16p . we have 

°= / ir + / 7r+ / 17 i= / 7r+/ 77 i + ° 1 - 2 - 18 
Jr 1 ov J r2 dv J Fa dv J Tl dv J T2 dv 

Meanwhile, by Green's formula and (|2.16[) . we have 



, ; dvi f dvi f dvt 

I Vliil = / V%— h / V\— h / V\- 



(x° 2 ,6)nn Jr 1 dv J T2 dv J T;i dv 

dvi f dvi f dvi 
dv ./ r „ dv J r . dv 



(2.19) 



Therefore, by combining ([2T7) . (|2~T8|) and (|2T9]l . we have 

012= / _ |V^| 2 +0(1). 

JB(x° 2 ,S)nri 

Similar to the energy estimates given in Lemma 1.5, Lemma 1.6, and Lemma 1.7 in [5], we have 

If C 

< / JVvi| 2 <^=, forn = 2 



Cy/ii2 jB(x° 2 ,S)nQ V^12 

hlne 12 \< [ JV V1 | 2 <C|lnei 2 |, for n = 3 

^ JB(x° 2 ,8)nn 

^ < [ _ |Vwi| 2 < C, for n> 4. 

JB(i» 2 ,j)nS] 
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Therefore, 



1 C 

< ai2 < , for n = 2, 



C v / £l2 i/£l2 

— | In £12 1 < 012 < C| In £12 1, for n = 3, 
— < 012 < C, for n > 4. 

□ 

Knowing enough properties of the system of linear equations (|2.10[) from Lemma 12.31 and Lemma 
, we have 

Proposition 2.2 Let € H (O) &e i/ie weaA; solution to equation and C, £/ie ua/we of on 

Di, then for any 1 < i 7^ j < to, i/iere exists a universal constant C > depending only on n, k q , r , 
119011(72,0, {||d_D,-||c | 2.<*}j independent of {eij} such that 



\d -Cj\ < C^/eJJ\\ip\\ L ^ {an) for n 

1 

In Si 



\Q -Cj\< Cji^— rlMU-tfn) /or n = 3, (2.20) 



|Ci - Cil < C\\ip\\ L oo {m) for n > 4. 

Proof: By Lemma T2. 31 we know that the matrix —A satisfies condition (Al) and (A2), then applying 
Proposition 1 2 . 1 1 on (|2.10j) . we have, for any 1 < i ^ j < m, 

C 

\Ci — Cj\ < — |MU°°(6«) 
a.ij 

where C is some constant depending on n, kq, ro, ||9f2||c 2 .°i {H^-DiHc 2 ^}? but independent of 

By Lemma l2.5[ we immediately finish the proof. □ 

Now we are ready to complete the proof of Theorem 11.11 
Proof of Theorem IJ.il We prove the estimates in dimension 2, the proof for the higher dimensional 
cases is similar. Without loss of generality, we assume i = 1, j = 2 and £12 < S. Now we need to prove 
the gradient estimates for Uqo in the narrow region between D\ and D 2 . For simplicity, we assume 

IMU°°(3fi) = !• 

By the decomposition formula (12. 7p . we have 



Vuoo = (Ci - C 2 )Vwi + C 2 (V(vi + v 2 )) + C iV v i + Vt, o 
By Lemma [2~4l we have 

C 

ll Vt, illL~(nnfl(aO «)) < — , ll V ^ol! L ^ ( onsK 2 , 5) ) < C ( 2 - 21 ) 

12 

where C is some universal constant. 



For i — 3, . . . , m, we have, by Lemma 12.41 

l|V«i|| i0 o ( n nB(a o aia)) < C. (2.22) 
Since v\ + v 2 — 1 on both dD\ and dD 2 , similar to the proof of Lemma |2~4"1 we can show that 

\\^(vi+v 2 )\\ L<xl(nnB{x o 2<5)) <C. (2.23) 
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By Proposition [Ml (E2Tj), (ET2"2"j) and (|2~2"3"j) , we have 

l|V« o|| i o 0( n nB(B o jiia)) < |Ci - ^HlVwill^^^^)) + |C a |||V(«i + «a)|| ioo(5nB(a o aia)) 

m 

+ X]l C,< lll Vu< lli~(nnB(x; 2 ,5)) + H Vu o|| L oo ( f 2ni 5 (;1 .o 2 5)) 

< C^iE— + C 
£12 

c 




As we mentioned in Remark 11.11 the strict convexity assumption of the two inclusions can be 
weakened. In fact, our proof of Theorem 1 1 . 1 1 applies . with minor modification, to more general inclusions 
as below. 

In W 1 , n > 2, for two closely spaced inclusions Di and Dj which are not necessarily strictly convex, 
assume dDi n B(0,r) and dDj n 5(0, r) can be represented by the graph of x\ — f(x') + ^j- and 
xi = -g{x') - e -f, then /(O') = fif(O') = 0, V(.g + /)(0') = 0. Assume further that 

\i\x'\ 21 <g(x') + f(x')<\ 2 \x'\ 21 , M\x'\<r/2, (2.24) 

where A 2 > Ai > 0, I G Z+. 

Under the above assumption, let E i? 1 (J7) be the solution to equation (|1.4[) . Then, for 
sufficiently small, we have 

_ ™=i 

l|V«oo||£o. ( nnB(x?^)) ^ ^IMU^an)^ 21 ifn-l< 21, 

llVttooll £ « M( onB( a »o f ,i)) < C lM^°°( 3n ) £ii | mE ..| ifn-l = 2l, ( 2 .25) 

l|Vuoo|| L oom nB(2: o 4)) < C||¥>||i« (a n}— i/n- 1 > 2i. 

where C is a constant depending on n, Ai, A2, ro, ||c?Z?i||c ,2 >° and ||9Z3j|jc72,a, but independent of ey. 
For the proof, please refer to the corresponding discussion after the proof of Theorem 0.1-0.2 in [5]. 



3 The insulated conductivity problem 

In this section, we consider the anisotropic insulated conductivity problem, which is described by 
Equation (|1.7p . As we mentioned in the introduction, the gradient can only blow up when two inclusions 
are close to each other. In order to establish the gradient estimates for this problem, we first consider 
the local version of the problem, namely Equation (|1.9j) . 

To make the problem easier, we first consider the equation in a strip. In this case, by using a 
"flipping" technique, we derive the gradient estimates in the strip. 
Denote, for any integer I 

1)5 < z x < (21 + 1)8, \z'\ < 1}, 
(2Z + 1)<5 and \z'\ < 1}, 
(21 - 1)5 and \z'\ < 1}, 

and 

Q = {zeR n ||z 1 | < 1 and \z'\ < 1}. 
We consider the following equation in Qo 

\d Zi (b i \z)d Zj w)=Q mQo, 
\ b 1] d Zj w = on T± 



Q, := {z eR n \(2l- 
T+ := {ze R n \z! = 
:= {z G W\zi = 
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where (b^) G C a (Qo)(0 < a < 1) is a symmetric matrix function in Qoj and there exist constants 
A 2 > A 2 > such that, for all £ G K", 

ll^'WIIc-Oo) ^ A 2- A 2|£| 2 < Vz eQ ,(e R". 

Then we have 

Lemma 3.1 Suppose w G -ff 1 (Qo) H £°°(Qo) a weafc solution of h3.1\) . then there exists a constant 
C > depending only on n, A 2 , A 2 , &ttt independent of 8, such that 

w/iere Qo(|) := {2 G E"||^i| < 6 and \z'\ < §}. 

Proof: For any integer 2, We construct a new function w by "flipping" w evenly in each Qi. We define 

w(z)=w((-l) l (z 1 -2lS),z'), VzeQ,. 

Therefore, we have defined w piecewisely in Q. 

We define the corresponding elliptic coefficients as follows 

for a = 2, 3, . . . , n, 

b al {z) = b la {z) = (-l) l b la {(-l)\z 1 -2l5),z l ), VzeQ,. 

for all other indices 

b ij (z) = b i '{{-l)\z 1 - 216), z'l Vz G Qi. 
Under the above definitions of w and b lJ , we can easily check that, for any integer I, 

d Zi (b ij {z)d Zj w)=Q in Q h 

b l3 d Zj w = onT±, 
Then for any test function ip G C^°(<2), we have 

= (by the definition of weak solution) 

Therefore w G i? 1 (Q) satisfies 

=0 inQ. 

Following exactly from |13j . we first introduce a new equation 

d Zi (B^(z) d Z] u) = inQ 

where 

{i im ze<3(. 2->((2i-i)i, o') & y (z) z £ Q h l>0; 
6 iJ (0) _ zeQo 

um zeQ<, *-»((2*+i)«, 0') &y ( z ) z e g;,i<0; 

then we define the norm 

||F||y.., - sup r 1 -'(-f r a\F\ p )' 

Q<r<l J 

Since b^{z) € C"(Qo) 1 b^ (z) is piecewise C a continuous in Q, then we can immediately check that 

\\b ij -B ij \\Yi+«,2 < C 
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where C is some constant only depending on A2. Using Proposition 4.1 in 13J, we have 

ll v ^lli-(iQ) < CNU 2 (Q) < C|Hk~(Q), 
Then by the definition of w, we have 

ll Vw lli-(Qo(i)) ^ C'lklU-CQo) 
where C > depends on n, A2, A2, but is independent of S. □ 

Since D\ and D2 are strictly convex domains, we can write 0(r), which is defined by (|1.8|) . as follows 

0{r) = {xeW l \- g(x') - e/2 < x x < f(x') + e/2, < r} 

With the side boundary r + and T_ as 

T+ = {x G R n |a:i = f{x')+e/2, \x'\ < r}, L={ie R n |a:i = -g(x')-e/2, \x'\ < r} 

where f(x') and g(x') are strictly convex functions, moreover they satisfy 

/(0') = 5 (0') - 0, V/(0') = Vg(O') = 0. 

Under the above notation, we prove Theorem 1 1.2 1 
Proof of Theorem ] Fix one point (0,x' o ) € C(§) and let 5 — y/ f{x' ) + g(x' Q ) + e, since f{x') and 
g{x') are strictly convex, then there exists a universal constant C depending only on ||9Z?i||c ,2 q and 
\\dD 2 \\c2.<* such that 

1 ; ' r' 12 



(3.3) 



-y\x Q r + s<6<C^\x' Q \i + e. (3.2) 

We shift the origin to (0, x' ) and rescale the coordinates with 8, then the new coordinates y = (yi, y') 
can be written as follows 

f Vi = xi/6, 
\y' = (x'-x' )/S. 

Let 

v(y) = u (6yi,x' Q + 5y'), a l \y) = a zj (Sy u x' Q + 8y'). 

Denote 

5(r) := {y G M"| - | - 5 K + V) < <% < f + /K + V), ly'l < *=} 
With its side boundary 

r+ := {y G R n |<$ yi = | + f(y' Q + 5y'), \y'\ < 7} 

f_ := {y G R n \S yi = -| - ff (y^ + V), |y'| < r}. 

By (13. 2p . we can find some universal constant r depending only on dD\ and 9Z?2, such that 0(r) is in 
the image of 0[r) under the above transform. Thus we have 

(d yi (a^d w v(v)) = in 0(r), 
\ a ij d y vv, =0 on f + U f _. 

where the coefficients a u satisfy, for some universal constant C, 

ll^'llc^f)) < C||a y ||c=(o W ) < CA U < ^ j (vMj (ty G 0(r), V£ G R"). 
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Next we construct a map $ : 0(r) i — > Qo, ®(y) = z with 



= , 6 yi +g{x' + Sy')+e/2 _ 
f(x' + Sy')+g(x' + Sy')+e ' 



(3.5) 



It can be verified directly that this map is a diffeomorphism from 0{r) to Qq. 
Let 

w(z) = v(^~ 1 (z)) 

Then from the definition of weak solution, we know that w(z) satisfies the following equation 

( d Zi (b^(z)d Zj w(z)) =0 inQo, 

\b l3 {z)d Z] w(z) = on r+ur . 



(3.6) 



where 



(&«(*)) 



|det(9^)| 

Therefore, we have transferred the original problem into Equation (|3.ip . 

In order to use Lemma ETTl we have to check that b l i(z) is strictly elliptic and || &* J " || c?™(Q ) ls t> oun d e d 
by some universal constant. First we show that there exists a universal constant A2 such that 



that 



A 2 |£| 2 V^eR", VzeQo 
Notice that the eigenvalues of (d y z) are A with multiplicity n— 1 and d yi z±. By 



1 2<5 2 



< C 



(3.7) 
we can prove 

(3.8) 



where C is some universal constant. 
Based on (13.81) .we have 



where A2 > is some universal constant 

The boundedness of ll^' 3 ll<7°(g ) can ^ e cnec ked similarly. 
Now applying Lemma |3. 11 we have 



II Vw IIl-(q (|)) ^ C \\ w \\l°°(Q ) 
Tracing back to uq through the transforms, we have, for any point x £ C(§), 

IV7 I M / C \\ u 0\\L^(O(r)) . C||uo|| L ~(0(r)) 
\VU (X)\ < < 



□ 



4 Appendix 

Some elementary results for the insulated conductivity problem 

Assume that in R", f2 and w are bounded open sets with C 2 ' a boundaries, < a < 1, satisfying, 
for some m < 00, 

m 
s=l 



16 



where are connected components of lo. Clearly lo s is open for all 1 < s < m. Given ip G C 2 (dil), 
the conductivity problem we consider is the following transmission problem with Dirichlet boundary 
condition: 

dxj{ (ka^(x) - a%(x))xu> + (x) d Xi u k ^ =0 in fl, 
Uk = <p on dQ, 



(4.1) 



where < k < 1, and Xu is the characteristic function of lo. 

The n x n matrixes A\ (x) := (a^(x)) in lo, ^2(2) :— {a^[x)) in 0\ZU are symmetric and 3 a constant 
A > A > such that 

A|C| 2 < a{{x)Uj < m 2 (Vi G w), A|C| 2 < a^z)^ < A|£| 2 (Vx G 

for all £ £ 1" and € C 2 (w), a 2 j (a;) G C 2 (H\w). 

Equation (|4.1[) can be rewritten in the following form to emphasize the transmission condition on 



d Xj (a\ 3 (x) d Xl u^j = 
d x . (a%(x) d Xi u^j = 



in lo, 
in f2\ZZJ, 

u k \+ = Uk\-, on duo, 

a£ {x)d Xi UkVj |, = fcay(x)9 Xi Ufci/j|_ on 



(4.2) 



[u k = p 



on. 



It is well known that equation (|4.ip has a unique solution in H (ft), and the solution Uk is in 
C 1 (f2\w) (~1 C 1 (aJ) and satisfies equation (|4.2[) . On the other hand, if Ufc G C 1 (f2\w) n C 1 (ZU) is a solution 
of equation (|4.2|) . then u k G i/ 1 (0) satisfies equation (|4.1|) . 

For G (0, 1), consider the energy functional 



h[v] 



a-l{x)d Xi vd Xj v + 



a? (x)d Xi vd Xj v, 



(4.3) 



defined on 



#i(f2) := {u G ^(0)1 w = <p on 90}. 



It is well known that for k G (0, 1), the solution u k of (|4.ip is the minimizer of the minimization 
problem: 

I k [u k ]= rma I k [v]. 



For fc = 0, the insulated conducting problem is: 

d Xj (a%(x) d Xz u ^j = 
a? (x)d Xi u Q Vj\ + = 

Uq = Lp 

d Xj (a^(x) d Xi uo) = 
t u |+ = ""ol-, 



in Q\lv, 

on 8lo, 
on dfl, 
in lo, 
on duj. 



(4.4) 



Equation (|4.4p has a unique solution uq G _ff 1 (f2), which can be solved in 17 \ w by the first three 
lines in (|4.4|) . and then, with uo|a w , be solved in u> using the fourth line in (|4.4|) . It is well known that 

uq g c 1 (f7 \ w) n c 1 (zu) . 

Define the energy functional 



IoM : = 



0\5 



a%(x)d Xi vd X] v, 



(4.5) 
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where v belongs to the set 

A Q := {v e fl" x (fi\5J)| v = tp on dfl}. 
It is well known that there is a unique vq £ Ao which is the minimizer to the minimization problem: 

I [v ) = mm I [v\. 

v£A 

Moreover, vo — uq a.e. in fi \ 57, where Uo is the solution of (|4.4[) . 
Now, we give the relationship between Uk and uq. 

Theorem 4.1 For < k < I, let Uk and uq in H 1 ^) be the solutions of equations \4->ty an d &4-4j) > 
respectively. Then 

Uk - 1 in H (f2), as k — > 0, (4-6) 



and, consequently, 

Proof: We will first show that 



lim Ik[uk] = Iq[uq]. (4.7) 



sup \\Vuk\\ L 2(n) < oo. (4.8) 

0<fc<l 



Since Ufe is the minimizer of in H^(fl) and «o := wo|o\u is t ne minimizer of Iq in .Ao, we have 
Afc k f ■ 

-yllVufclU 2 ^) + < q / tt i J W 9 i,' 1 t 9 i J ,| i: + - f oW 

fc f ■ 

< - / ai J (aj)9 a;i Ufc9 a . i Ufe + /o[ufc|n\s?] = 4M 

< /fc[«o] = 2 / a i ( x )9xiU d X:j uo + J N]) 
Afc 

< -g-HVitoll^Co)) + Jo[«o]- 



Thus 



SUp ||VUfe||L2( w) < oo. 
0<fc<l 



On the other hand, 



^||Vw fc ||L 2 (o\^) < h[uk] < h[uo] < yII Vu oII£ 2 (o)- 
Estimate (|4.8[) follows from the above. 

Since Uk — <p on 9f2, we derive from (|4.8p that sup 0<fc<1 ||ufc||iji(o) < oo. Let Uk — v wjj fe Tf^fi) 
along a subsequence of fc — > (still denoted as fc — » 0). 

We will show that Uq is a solution of equation (14. 4ft . Therefore, Uq = uo- 
We only need to establish the following three properties: 

d Xj (ji(x) d x .u^j = in n\57, (4.9) 
^ (a? (a?) =0 few, (4.10) 

Uq € C 1 (ri \ uj), a l 2 {x)d Xi UQVj I = onduo. (4-11) 

(i) For fc € (0, 1), we see from equation (|4.1I) that 

d Xj (a% (x) d Xt u k ^ = 0, in Q \ UJ, 
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d Xj laf (x) d Xi Ukj = 0, in lu. 

Since u k converges to Uq weakly in (|4. 9[) and (|4.10[) follow from the above. 

(ii) For any w £ Aq, we extend it to w G H^(Q) (i.e. w = w in fl\uJ). By the minimality of Uk, 

Ik{uk) < h{w). 

Sending A: to leads to 

Io{v>l\n\u>) < Iq{w). 

Thus Uq = u a.e. in O \ u. (|4.11|) follows. 

We have proved (14. 6p . Theorem 14. II is established. 
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